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SUMMARY 


The success of future permanent space stations depends on the 
development of a space shuttle vehicle having aerodynamic maneuvering 
capability. The purpose of this technical report is to investigate the 
optimum maneuver of such a vehicle reentering a spherical, stationary, 
and locally exponential atmosphere. The use of Chapman's modified 
variables and a rescaled lift-drag polar leads to the formulation of a 
set of dimensionless equations of motion for flight analysis. The 
resulting equations are exact in the sense that they are also valid for 
flight in the vacuum. For the vehicle, we only have to specify the most 
important performance parameter, namely the maximum lift-to-drag 
ratio E . On the other hand, the planetary atmosphere is characterized 
simply by the so-called Chapman's atmospheric parameter pr. 

For planar flight several typical optimum maneuvers are investi- 
gated at different altitude ranges, low, moderate and very high. In each 
case the characteristics of the optimum lift control are discussed. 

For three-dimensional flight, the procedure to solve the optimum 

trajectory for maximum cross range is discussed in detail. Finally, 

using the equilibrium glide condition the maximum cross ranges for 

^ 

entry from circular speed, for several values of E , and the footprint 
for E* =1.5 are computed in this reduced problem. A technique of 
coordinates rotation is used which makes the iteration procedure for 
solving the footprint of a reentry vehicle much more effective and 
geometrically meaningful. 


XI 



CHAPTER 1 


INTRODUCTION 

The purpose of this report is to investigate the optimum 
maneuver of a space vehicle having aerodynamic maneuvering capa- 
bility. The present day space shuttle is an example of such a kind of 
space vehicle. In the previously published literature, the analyses 
are either for constant lift-to-drag ratio [l,?], or usually numeri- 
cally oriented and is confined to the performance of a particular 
vehicle [s] . In order to maintain the generality of the results, we 
shall introduce a set of dimensionless variables and a rescaled lift- 
drag polar, to derive the dimensionless equations of motion for flight 
analysis inside a spherical, stationary, and locally exponential 
planetary atmosphere. The resulting equations are exact in the sense 
that they are also valid for flight in a vacuum and are almost free 
from all the physical quantities of the vehicle and the planetary 
atmosphere. For flight at very high altitude with orbital speed, a 
Newtonian, inverse- squared force field is used. By a simple canon- 
ical transformation [5,6], the corresponding eqriations for low 
altitude and low speed flight over a flat earth model are obtained. 

Two naain types of optimum maneuver in a vertical plane vdll 
be investigated at three different altitudes, low, moderate, and very 
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high. In the pull-up type maneuver we either maximize the final speed 
with the final altitude prescribed, or vice versa. At very high 
altitude with orbital speed, the maneuver generates the useful skip 
trajectory. In the gliding type maneuver we maximize the gliding 
range. The three-dimensional gliding maneuver for maximum cross 
range will also be discussed and then solved in a reduced problem. 
The footprint of a reentry vehicle will be assessed. 

The organization of the report is as follows. After this intro- 
ductory chapter, the dimensionless equations of motion for three- 
dimensional atmospheric flight are derived in Chapter 2. The prob- 
lem is then formulated as an optimal control problem with the adjoint 
equations and the control law derived in Chapter 3. In Chapter 4, the 
equations for planar flight are deduced from the general equations of 
Chapter 3. Then in Chapter 5 they are transformed into the eqtiiva- 
lent form appropriate to the flat earth model. The numerical applica- 
tions are carried out in Chapters 6, 7, and 8. The case of flight over 
a flat earth is first analyzed in Chapter 6. The concept of a linear- 
ized singular arc [s] is introduced and tested. The case of planar 
flight over a spherical earth is discussed in detail in Chapter 7. Then 
in Chapter 8 we discuss the procedure to solve the three-dimensional 
optimum trajectory. The problem is then simplified and solved with a 
footprint obtained. The final chapter, Chapter 9, summarizes the 


main results obtained. 



CHAPTER 2 


DIMENSIONLESS EQUATIONS OF MOTION 

In this chapter, the three-dimensional equations of motion of 
a nonthrusting, lifting vehicle entering a stationary spherical plane- 
tary atmosphere are introduced. Then by using the modified Chap- 
man's variables, a normalized lift coefficient, and a dimensionless 
arc length as the independent variable, a set of dimensionless state 
equations are obtained for entry analysis. It will be seen that, by 
this formulation, the only physical parameter involved is the maxi- 
mum lift-to-drag ratio, and the planetary atmosphere is simply 
characterized by a value referred to as Chapman's atmospheric 
parameter. 


2, 1 Three-Dimensional Equations of Motion 

The equations of motion of a nonthrusting, lifting vehicle 


entering a stationary spherical planetary atmosphere are 

^ = V sinV 
dt 

dV pSCoV^ 

„ dY pSClV^ I V^\ 

^ ^ (s - 7 "; 


3 




4 


d9 V cos Y cos 
dt r cos 4> 

d4> _ V cos Y sin 4^ 
dt r 

d4i rSClV^ . y2 

V — = rr sin <r - — cos Y cos 4J tan <b 

dt 2m cos Y r 

where t is time, (r, V, Y, 9, cj>, 4*) are state variables and are defined 
in Figure 1, p is density of the atmosphere, S is the reference area 
of the vehicle, Cj-j and Cj_^ are the drag and lift coefficients, m is the 
mass of the vehicle, g is the magnitude of gravitational acceleration, 
and 0 - is the bank angle. The flight path angle Y is defined to be 



Figure 1. State Variables, Control Variables, and Other 
Parameters Defined with Respect to Inertial 
Coordinate, OXYZ. 
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positive when the velocity is directed above the horizontal plane. The 
bank angle tris taken to be positive for a bank to the left. For flight 
in a Newtonian force field of a spherical planet, the magnitude of the 
gravitational acceleration is of the form 




( 2 . 2 ) 


where p is the gravitational constant. The density of the atmosphere 
p is assumed to be locally exponential, that is, it obeys the differen- 
tial law 

A n 

(2.3) 


= - p dr 


where the inverse scale height p is a function of the distance from the 
center of the planet r. 

There are two control variables, one is the bank angle or and 
the other is either the lift coefficient or the drag coefficient Cjj. 
For a given vehicle there is a lift-drag relation; therefore either the 
lift coefficient or the drag coefficient Cq can be used as the con- 
trol, We shall use as lift control a normalized lift coefficient \ such 


that 

(2.4) 


where is the lift coefficient corresponding to maximum lift-to- 
drag ratio E"". If is the corresponding drag coefficient, then 

Cj^ = C^U\) (2.5) 

where f(\) is the ftmction specifying the chosen drag polar (see 
Appendix A). When X = 1, the flight is at maxim\im lift-to-drag ratio. 
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Thus we also have f(l) = 1. We shall consider a parabolic drag polar 
with the simple function 

f{\) = i(l+X^) (2,6) 

In general, the parameters Cj^, , and E'‘' are functions of Mach 
number; but in the hypervelocity regime they are essentially constant. 


2.2 Dimensionless Equations of Motion 

The following dimensionless variables are introduced. 


Z = 


PffL 

2m 




V = 


gr fi/r 
t 


(2.7) 


/ V 

— cos y dt 

r 

0 

where Z and v are the modified Chapman's variables [l], Z is pro- 
portional to the atmospheric density p and will replace the altitude, 
while the dimensionless kinetic energy v is a measure of the speed. 
The remaining dimensionless variable s is the dimensionless arc 
length. It is monotonically increasing and will replace the time as 
an independent variable. By using Eqs. (2.4), (2,5), (2. 6) and (2, 7) 
in Eqs. (2. 1), we have the dimensionless three-dimensional equations 
of motion 

^ = - k^Z tanV 
ds 

dv _ k Z V (1 + X^) 
ds E* cos V 


- (2 - v) tan Y 
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ds 

ii 

ds 
ds 


k Z X cos cr 
cos Y 




cos 


4^ 


cos <j) 


= sin 4 j 


( 2 . 8 ) 


d^* k Z X sino- . . , 

T“ = ft; — - cos ^ tan cj) 

ds cos Y 

where E is the mcLximum lift-to-drag ratio, and k is the dimension- 
less product §r. These equations are exact and hence are valid for 
Keplerian motion outside the planetary atmosphere. The only slight 
simplification is that in the equation for Z, the exact coefficient of 
- Z tan Y is 

1 < 


k^ = k^ 


(- 


2k^ Zp^ dr / 


( 2 . 9 ) 


For a strictly exponential atmosphere, p = constant and 

^=0 
dr 

On the other hand, if an isothermal atmosphere is considered, 
P/g = constant and 

-L_ ig. = -L 

Zp^ dr ■ k^ 


( 2 . 10 ) 


( 2 , 11 ) 


t-2 


In both cases k~ is a function of k^ = (3r. Chapman has shown that in 
the reentry range of the altitude, this product is oscillating about and 
near a mean value [l]. Furthermore, its value is much greater than ' 
unity, e. g. , for the earth's atmosphere k^=i 900, thus we take 


k^ k^ 


( 2 . 12 ) 
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The Eqs, (2,8) are the state equations for entry analysis. It 
is seen that the only physical parameter of the vehicle involved is the 
maximum lift-to-drag ratio E*. Furthermore, any planetary atmo- 
sphere is simply characterized by a properly selected value k^. This 
mean value will be referred to as Chapman's atmospheric parameter. 



CHAPTER 3 


VARIATIONAL FORMULATION 


With an adjoint vector introduced, we formulate the problem 
as an optimal control problem by using the Pontryagin's maximum 
principle. The control law is derived. The integrals of the motion 
are obtained. Then there is a change in the adjoint variables to have 
a better form for the adjoint equations. Finally, the parameters of 
the problem under different cases are discussed. 


3. 1 Variational Formulation 

The Eqs, (2,8) are the state equations with two control 
variables, the lift control X and the bank angle cr. They are subjected 
to the constraints 


|\1 <X 

|(T I < 0- 


max 


(3, 1) 


' ' max 

These controls are to be selected to bring the vehicle from a certain 
prescribed initial condition to a certain partially prescribed final 
condition, such that a certain function of the final state variables is 
minimized. 

Using the maximum principle, we introduce the adjoint vector 
p to form the Hamiltonian 


9 


10 


H = - k= Zp^tanY - P, 

fkZXcoscr . 1 \1 . cos 4^ , . 

+ p + (1 - — ) + T + P , sin 4 j 

t L cos y X V / J 


P 4 . 


0 cos cj) ^ 4 > 
"k Z X sine 


cos^ V 


cos 4* tan c|) (3.2) 


where p^, x = Z , v, Y, 0 , <|), and <4 are the adjoint components corre- 
sponding to the six state variables, respectively. They are governed 
by the following adjoint equations 


^Pz 

ds 


,7 . k v(l +X^) kX cos e 

= k p ^ tan Y + p -rr* ^ - p ~ 

^Z ^v E*cosY ^Y cosY 


k X sin cr 

5 I 2 

4 * cos *Y 


ds 


I E* 


( 1 +X^) 
cos Y 


tan Y 


]- 


Zl 

v" 


dp. 


ds cos^Y 


,7„ rkZv(l + X^) sin Y 

k^Zp +p ^ 

Z V L 


E* 


+ (2 


-V,] 


-kp ZXcos esinY - 2kp, ZXsine tanY 
Y 4* 


Zfe 

ds 


= 0 


(3.3) 


dp 


i 


ds 


cos 4 r . , 1 

TT L-Pa sin 4 + p, J 


cos 4 ^9 


4^ 


dP4 

ds 


sin 4 


cos 


5 tpg - P 4 Sin4] - P^ cos 4 


The solution is then obtained by integrating the two sets of state and 
adjoint equations, subjected to the end conditions, and at each instant 
selecting the lift control X and the bank angle cr such that the Hamil- 
tonian is an absolute maximum. 
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The Hamiltonian H will be maximum either at the boundary of 
the control set or at an interior, variable point where. 


9H 

8 X 


= 0 


9H 
9 cr 


= 0 


(3,4) 


Explicitly, we have 


X cos (T = 


Py 


2 vp 


\ sin 0 * = 


TT* 

E P4; 


2 vp cos Y 

V 


(3,5) 


3, 2 Integrals of the Motion 

It is known that the problem has a number of integrals [2-4] . 
First of all. 


H = Co (3,6) 

where Cq is a constant. Then by solving the last three equations of 
Eqs, (3,3), we have 


^0^ 


= Cl 

= C^cosQ + C 3 sin 9 

= Cl sin cj) - cost}) (C^sinS - C 3 COS 0 ) 


(3,7) 


where Ci , C^ , and C 3 are constants of integration. 

To simplify the first three equations which are not integrable 
analytically in Eqs. (3. 3), it is convenient to use the modified adjoint 
variables defined as 


P = k^Zp^ 


N = vp^ 


Q = p 


(3.8) 
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The corresponding modified adjoint equations are then 

_ 2 


dP _ k^Z r N 
ds cos "Y Le* 


(O' 




cos^ y 


)] 


^ - (Q + ZNtanY) 

ds V 


(3.9) 


dQ ^ 1 

ds cos^Y 


kZsinY (2 - v) ] E*k Z sinY f 2 


p + nl— p— + 


4N 


L cos Yj j 


In terms of P, N, and Q, the optimum lift and bank controls become 


X cos a = 


X sin cr = 


E*Q 

2N 

E* 




2Ncos Y 


2 E*^ / 2 ^4^ 

X^ = — 


4N^ 


cos^ Y, 


and the Hamiltonian becomes 


(3.10) 


[5^ 


■PtanV - N I ^ tanvl - 

cosY V J V 


+ 


E"^k Z 
4N cos Y 


(q^ 


Jj 


+ P, 


cos 


4^ 


+ p , sin 


cos“^Y^ ^9 cos (j) ^ 4 > 

- p cos 4 ^ tan 9 = Cq 
9 


(3. ii; 


In summary, the optimal solutions of this problem are 
governed by the Eqs. (2.8) for the state variables, Eqs. (3, 7) and 
(3.9) for the adjoint variables, and Eqs. (3. 10) for the controls. It 
requires six parameters Ci , C 2 , C 3 , P^, N^, and to satisfy the 
final and transver sality conditions, where P. , N. and Q. are the 
initial values of P, N, and Q, respectively. The Hamiltonian equation, 
Eq. (3.11), can be used to check the accuracy of the integration. 
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For the most practical cases, the arc length s is not 
prescribed at the final time. Thus H = Cg = 0 and Eq. (3. 11) becomes 


P tan V - N 


k Z 


E* cos V 
E*k Z 


V 


(1 - V) Q 

V 


4Nc 






cos^ Y 


\ cos 4^ . I 

j + Pa TZ7T + Pa sin 4^ 


0 cos 4> 


4j 


- p cos 4i tancj) s 0 

4 ^ 


(3.12) 


Using this integral, one of the three equations in Eqs. (3. 9) can be 
deleted. But there are some difficulties in so doing. First, Eq. 

(3. 12) is quadratic in both Q and N. To solve either Q or N from 
Eq. (3, 12) requires frequent change in sign in front of the square 
root each time the quantity iinder the square root passes through the 
value zero. Next, in Eq. (3, 12) the coefficient of P is tanY. When- 
ever Y goes to zero, P cannot be determined. Hence, it is more con- 
venient to use Eq. (3. 12) solely to determine one of the three initial 
values, either P. or N. or Q. , and to check the accuracy of the inte- 
gration. Anyway, it is obvious that the number of parameters is 
reduced by one, that is, from six to five. 



CHAPTER 4 


PLANAR FLIGHT 

In this chapter, we deduce the governing equations for the 
optimum reentry trajectories confined to the plane of a great circle 
from the general three-dimensional equations of Chapter 3. They are 
the state equations, the adjoint equations , the control law, and the 
Hamiltonian integral. Then by a change of adjoint variables, we 
obtain a handy equation for the control variable \, and the number of 
parameters is reduced by one. 


4. 1 Governing Equations 

For entry trajectories in the plane of a great circle, we have 

0- = <j> = = 0 (4, 1 ) 

and the independent variable s is simply the range angle 9. The state 
equations and the modified adjoint equations are reduced to 


dZ 

d9 


Z tan y 


dv _ k Z V (1 -f X^) 
d9 E"^cosV 


(2 - v) tan Y 


^ kZ\ ^ N 

d9 cos Y V V / 


(4.2) 
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and 

dP _ k^Z r^L E*Q^ 1 
d0 cos y L E* J 


= - - (Q + 2N tanV) (4.3) 

d0 V 


dQ 1 

jp + N 

'k Z sin Y ^ (2 - v) ' 

E*kZ Q^sinY 

d9 cos^Y 

E* V 

4N 


respectivelyo 


The optimum, lift control is either \ = or a 


variable X ,such that 


E* Q 


(4. 4) 


The Hamiltonian integral becomes 


^PtanY - N 


kZ ^ (2-v) ^ ^ 

+ ^ tanY 

V 


E*cos Y 


J- 


(l-v)Q + E 'k Z Q" ^ 

V 4 N cos Y ^ 


where Ci is the same constant of integration as in Eqs. (3,7). 

In general, this is a three-parameter problem, with P^, N-, 
and Qj as the three parameters. For the special case where the 
range angle 0 is not prescribed at the final time, i. e, , 0£ is free. 

Cl = 0 and it becomes a two-parameter problem. 


4, 2 Change of Adjoint Variables 

From the expression of Eq, (4.4), it is seen that a simplification 
can be made if we use X as a new variable. Another variable which 
will be seen to be useful is 



(4,6) 


Using (F, N, X) as a new set of variables to replace the modified 
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adjoint variables (P, N, Q), and taking the derivative of Eqs, (4.4) 
and (4. 6) with respect to 0, we have 


dF _ Z (1 -X^) 2F 
d0 “ E*cos V ^ E* V 


(X + E'” tan Y) 


=-^(X + E*tanY) (4,7) 

dX _ kZ (1 -X^) sinY 2X (X + E*tanY) E* / 
d0 2cos^Y ^ E* V 2cos^Y \ v/ 


The Hamiltonian integral Eq. (4. 5) in terms of (F, N, X) becomes 


kZ (1 -X^) 2 (1 - v)X 

E*cosY ^ E*v 


+ 




tan Y 



(4,8) 


In Eqs, (4,7), it is seen that the first and third equations are inde- 
pendent of N, It can be shown that N^, the initial value of N, is free 
whenever the final value of N doesn't appear in the trans versality 
conditions. Thus the second equation of Eqs, (4,7) can be deleted. 

It becomes a two-parameter problem for the general case. For the 
special case if 0^ is free, Cx = 0 and it is simply a one-parameter 
problem. The Hamiltonian integral for this special case is, from 


Eq. (4.8), 


k Z (1 -X^) 
E* cos Y 


E* V 


+ 




tanY = 


0 


(4,9) 


As has been mentioned in Section 3,2, there are difficulties in 
using Eq. (4, 9) to solve for X or F, To solve for X from Eq, (4. 9), 
we have to determine the sign in front of the square root and change 
this sign each time the quantity under the square root passes through 
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zero. At that instant the equation has a double root 

(1 - v) cos Y 


(4.10) 


From the third equation of Eqs. (4, 2), it is seen that this corresponds 
to dY/d0 = 0. Physically, the flight path angle passes through a maxi- 
mum or a minimum and the trajectory has an inflection point at this 
instant. This behavior is typical in an optimal trajectory. Therefore 
it is more convenient to obtain the optimum X directly from integra- 
tion. On the other hand to solve for F from Eq, (4. 9) will become 
impractical whenever Y is approaching and passing through the value 
zero. Hence, Eq. (4, 9) will be used solely to compute the initial 
value Fj in terms of X^ and to check the accuracy of the integration. 



CHAPTER 5 


FLAT PLANET SIMPLIFICATION 

The equations we have derived in the preceding chapter are 
the optimum equations for the general case of planar flight. They are 
to be used when the speed of the vehicle is of the order of orbital 
speed, v 1 , which occurs at high altitude where the value of Z is 
small. They are, of course, also valid at low altitude and low speed. 
But in this case, without compromising the accuracy, it is simpler to 
use the equations within the framework of a flat planet model. These 
equations are to be deduced in this chapter. 


5, 1 Governing Equations for Flat Planet Model 

It is interesting to know that by a proper change of variables 
we can deduce the dimensionless equations for the flat planet case 
from the general equations of planar flight in the preceding chapter. 
At low speed and low altitude, it is more convenient to use the 
following dimensionless variables 


w 


2 m P 
pSC^* 



PY 


(5, 1) 


where w is the dimensionless wing loading which will replace the 
altitude, u is the new dimensionless kinetic energy to represent the 
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speed, and y is the dimensionless linear downrange. The relation- 
ships between the two sets of dimensionless variables (Z, v, 0) in 
Eqs. (2,7) (where s has been replaced by 9 in planar flight case) and 
(w, u, y) in Eqs. (5. 1) are 


Z 

w 


_u 

k' 


9 - 

0 -k2 


(5.2) 


Since the value of k^ is much larger than u, e. g. , for the earth's 
atmosphere k^ ^ 900 and u is of the order of unity at low speed, we 
have 


k^ » u 


(5.3) 


This is the flat planet condition. Upon substituting Eqs. (5. 2) into 
Eqs, (4.2) and using Eq, (5.3), we have 

dw 


dy 

du 

dy 


= w tanV 


u(l +X^) 
E*'' w cos Y 


- 2 tanY 


(5.4) 


3 X 1_ 

dy wcosY u 


These are the state equations for flat planet model. We will obtain 
identical equations by starting out from the classical equations for 
flight over a flat planet and using Eqs. (5. 1) in them. It is seen that, 
although an exponential atmosphere is still used for this case, the 
characteristic parameter k^ of the atmosphere is removed from the 
equations. Hence, the flight behavior is independent of any particular 
atmosphere. 
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Again, we can use Eqs. (5.4) to form the Hamiltonian and 
derive the optimum equations, as has been done in Chapter 3. It is 
more elegant and informative to use the condition of a canonical trans- 
formation as a handy tool to effect the transformation from the old to 
the new variables [5,6], This, coupled with the condition of Eq, (5.3), 
will lead directly to the equations for the optimal control of the flat 
planet case. 

For a transformation from the variables (Z, v, 9) with 
Hamiltonian H to the new variables (w, u, y) with the Hamiltonian 
to be canonical, we have the necessary and sufficient condition that 
the quantity 

(p„ dZ + p dv - H d0) - (p dw + p du -«3^dy) = dU (5. 5) 

Zu V "W VI 


be an exact differential. In particular, for dU = 0, and using Eqs, 
(5. 2), we have 

H = 


w^Pw 

Pz = -^ (5.6) 

Pv = ‘'"Pu 

p,^ remains \mchanged. By using Eqs, (5, 2) and (5. 6) and the con- 
dition (5.3) in the Hamiltonian integral (4, 8), we get the Hamiltonian 
integral for the flat planet case 

1 - 2X /_ 2 


+ - (g - -) tanY = — 

E u V u^ u 


(5.7) 


u 


E"'' w cos Y 

where again Cj = is a constant of integration with Ci = 0 for the 
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free range case. The variable G in Eq. (5.7) is the analogue of the 
variable F in Eq. (4. 8), and is defined as 


wp. 


G = 


w 




(5.8) 


Similarly, performing the same transformation and using the same 
condition on the first and third equations of Eqs. (4.7), we have 


dG 

dy 

6X 

dy 


+ E*tanV) 


E^ w cos y 

(1 - X^) tanV 
2 w cos Y 


+ 


E^u 

2 \ 

E • u 


(X + E'-'tanY) 


2 cos^ Y 


(°-f) 


(5.9) 


Again, this is a two -parameter problem in general. It will be 
reduced to a one-parameter problem when y^ is free and Cj = 0. The 
Hamiltonian integral for this special case is 


1-X' 


2X 


E^wcosY E* 


- - (g - -) 

U ^ U / 


tanY = 0 


(5.10) 



CHAPTER 6 


OPTIMAL TRAJECTORIES FOR FLAT EARTH 

In this chapter, two categories of optimum trajectories are 
computed numerically using the equations derived in the preceding 
chapter for a flat planet model. The first category of optimum tra- 
jectories is for the pull-up maneuver. We either maximize the final 
speed with a prescribed final altitude or vice versa. The final flight 
path angle can be either prescribed or free. We consider both cases 
of unconstrained X and constrained X. The second category is for the 
glide trajectory which maximizes the final range with prescribed final 
altitude, final speed, and/or final flight path angle. Since the equations 
used are independent of the planet and its atmosphere, so are the 
results. But to have some idea about the physical quantities of the 
flight, we use the flat earth model and its atmosphere as an example 
to get dimensional q'uantities from the dimensionless results. In the 
last section of this chapter, the linearized singular control technique 
is introduced and tested. It is proved to be useful in the saving of the 
computational work. 
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6. 1 Maximum Final Speed or Maximum Final Altitude 

In this case, it is proposed to find the optimum lift control to 
bring the vehicle from the initial condition 


y = = 0 , w = , u = u^ , "Y = (6, 1) 

to the condition at the final instant y^ such that either 

w = , u = u^ = maximum (6, 2) 

or, u=u^,w = = maximum (6.3) 


We call this the pull-up type maneuver. A sketch of this type of 
trajectory is presented in Figure 2. The condition of Eq. (6.2) is to 
maximize the final speed with a prescribed final altitude, while the 



Figure 2. Geometry of a Pull-Up Maneuver 
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condition of Eq. (6, 3) is to maximize the final altitude with a pre- 
scribed final speed. They are eqtii valent, and it will be shown that 
their solutions are obtained through a single formulation. We shall 
assume that the final range is free, thus Cj = 0, and thus we can use 
Eq, (5, 10). It is a one-parameter problem. 

As has been explained before, to avoid the difficulties in using 
Eq. (5.10) to solve for either X or G, we shall integrate both of the Eqs, 
(5,9) along with the state equations, Eqs. (5.4). For the five initial 
values required, since the initial state (w^, u^, Y^) is given, we need 
only the two initial values and to start the integration. We set 
to be the only parameter of this problem, and obtain G^ from Eq. 
(5. 10). This can be done except when Yj = 0. The case with Y^ = 0 
will be discussed later in this section. 

For the numerical computation, we shall use the initial state 

(w , u. , Y.) = (, 5, . 5, - :^) (6. 4) 

Although a specific set of values has been used, it is found that the 
optimum lift control has a general typical behavior. For the maxi- 
mum lift-to-drag ratio E^‘', we shall use E"'' = 10 which is typical for 
a fighter aircraft, and E"" = 4. 5 which is somewhat higher than the 
value of a shuttle vehicle at low speed. To maximize the final speed 
with a prescribed final altitude, we start the integration with a 
guessed and stop it at w = w^. If the final flight path angle Y£ is 
prescribed, this value is used to adjust X^ xintil the condition is met. 



25 


The resulting trajectory is the optimal for majximizing u^ with the 
prescribed W£ and satisfied. If is free, then by the transvers- 
ality condition, p = 0. From the first equation of Eqs. (3, 5), after 
being transformed to the form for the flat planet case, we have at the 
final instant 


E*p, 


Xr - 




= 0 


(6.5) 


This condition is used to adjust for the free case. The result is 
the overall best since the final flight path angle is also optimized. 

A similar procedure is used to find the optimum trajectory for 
the case of maximum final altitude with a prescribed final speed. 

Since the problem has one arbitrary parameter, namely the 
initial value X., the family of optimum trajectories is generated by 
simply integrating the Eqs. (5.4) and (5. 9) for different values of X. 


until X^ = 0. The results for E *' = 10 and 4, 5 are presented in Figure 
3, which is plotted in the ratio w/w. versus the ratio V/V^. The 
solid lines are the different optimum trajectories leading to the 
terminal boundary represented by the dashed line. From the defini- 
tion of w in Eq, (5, 1), if an exponential atmosphere is used, the 
actual altitude change is simply 


1 /^f 

Ah = h - h. = - log 1 — 

f 1 p ®\w.^ 


( 6 . 6 ) 


For any prescribed change in altitude, we can evaluate the cor re-' 
spending minimum speed reduction along the dashed line. Conversely, 



ure 3. Solution 
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we can evaluate the maximum altitude gain if w^/w^ > 1 , or the 
minimum altitude loss if w^/w^ < 1, for any prescribed speed 
reduction. 

Although the figure is plotted for a specific initial state given 
in Eq. (6,4), the use of dimensionless variables allows a general 
discussion of the influence of different physical characteristics of the 
vehicle on its performance. For a numerical example, with u^ = . 5 
and taking g = 9.81 m/sec^, l/j3 = 7162 m, the initial speed is 187.43 
m/sec or 674.7 km/hr. Assume a prescribed reduction in the speed, 
say = .7. Then flying optimally, the maximum final altitude 

is identified in Figure 3 along the dashed line of E* = 10 to be 
Wf/Wf = 1. 07. From Eq, (6. 6) this represents an altitude gain of 
484. 57 meters. The initial altitude with w^ = . 5 is 



which is a function of the wing loading m/SC *. For a higher wing load- 
ing, the same gain in the altitude can only be achieved at a lower alti- 
tude. In other words, small wing loading favors the pull-up maneuver. 

Figure 4 presents the variation of the normalized lift coefficient 
\ as a function of the flight path angle "Y for several optimal trajec- 
tories, Higher values of correspond to smaller speed reductions. 

It is interesting to notice that when 7 = 0, that is, when the vehicle is 
at the lowest point (or bottom) of the trajectory, the \ for different 
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Figure 4. Variation of the Optimal Lift Coefficient in Terms of the Flight 
Path Angle in Pull-Up Maneuver 
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trajectories with the same E'^ have nearly the same value From 

the Hamiltonian integral Eq. (5. 10) with Y = 0, we have 


2 \. 


u 


w/b >^b " ^ 


( 6 , 8 ) 


On the other hand, from the definition of u and w in Eqs. (5, 1), we 
have 

n 1 / « \ / S C T* \ 

( 6 . 9 ) 


_ 1 f P V 

W “g \ 2 


m 


Therefore, if ^-b i^ nearly the same for all trajectories, the corre- 
sponding dynamic pressure (1 p V^)^ is nearly the same, which in turn 
means that the indicated speed at the lowest point is nearly the same. 
Furthermore, at the lowest point, the normal acceleration as felt by 
the pilot is the opposite of the acceleration due to the lift force, which 


in terms of is 






2X^ 
b b 


1 -'b ‘ . 1 


( 6 . 10 ) 


Thus it is also nearly the same for all trajectories regardless of the 
final condition achieved. 

It is possible to obtain an approximate analytical expression 
for X^ by considering a particular trajectory in Figure 4 which shows 
a near constant value of X from Y = Y^ to Y = 0. From the second 
equation of Eqs. (5. 9), since dX/dy 0 at Y = 0, we have 

E* (uG - 2)^ = (6.U) 


Secondly, from the Hamiltonian integral Eq, (5. 10) at the initial 
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instant, 

VL^(l.-Xl) 

(uG-2), tanV. = r- + 2\. (6,12) 

1 1 w. cos Y. 1 

1 1 

It is confirmed by the numerical results that the product uG also 
varies slowly. Thus 


X. , (uG- 2). (uG - 2X (6.13) 

lb 1 b 

Combining Eqs. (6. 11), (6. 12), and (6. 13) gives the approximate 
equation for evaluating X^ 


_(w)j + ^ 


X^ - 2 X, cos Y 




(6.14) 


The values of X^ obtained from this equation is in excellent agreement 
with the numerical results, as shown in Table 1. 


Table 1. Comparison of Approximate and Actual X^ 
(u/w)^ variable, = - 1/E* 



(u/w). 

1 

. 6 

. 8 

1.0 

1.2 

1.4 

E''^ = 10 

Approx. X^ 
Actual X^ 

3. 833 
3. 825 

2.972 
2. 964 

2.491 

2.487 

2. 188 

2. 186 

1. 982 
1.981 



(u/w)^ 

. 6 

.8 

1.0 

1.2 

HI 

E* = 4, 5 

Approx. X^ 
Actual X^ 


3. 598 
3. 635 

2. 861 
2. 879 

2, 434 
2. 456 

2.159 
2. 172 
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To assess the influence of the initial speed,- we use the same 
values = . 5 and = -1/E , and generate several families of opti- 
mum trajectories using u^ = .3, .4, ,5, .6, and ,7. The solutions, 
that is, the terminal boundaries of different families, are presented 
in Figure 5. It is obvious that higher altitude gain is obtained with 


higher initial speed. 
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Finally, using the same values = , 5 and = .5, we vary, 
the Yj to analyze its effect on performance. The solutions are pre- 
sented in Figure 6, Obviously, the performance improves as Y. 
increases and becomes positive. One interesting observation is that 
when Yj^ = 0, that is when the maneuver starts horizontally, can be 
solved from Eq. (5. 10), 
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Thus it is solely a function of the initial dynamic pressure (u/w)^ 

(or the initial indicated speed), and has the same value for all tra- 
jectories. In the example here its value is 2.4142. This also means 
that can no longer be a parameter in this case. To generate the 
family of optimum trajectories, we have to use either or (dX/dy). 
as a parameter. 

6. 2 Pull-Up Maneuver with Bounded X 

In the preceding section, to display the behavior of the lift 
coefficient along an optimum trajectory, we put no restriction on its 
upper limit. This is of no problem for vehicles with high maximum lift- 
to-drag ratio since the optimum X is within a reasonable limit. But 
for vehicles with low maximum lift-to-drag ratio, as in the case of 
the reentry vehicle, the optimum X may be unacceptable since it can 
exceed the stalling lift coefficient X . 

rxj.AX 

To discuss the behavior of the optimum trajectory in the case 
of bounded X, we refer to Figure 7 which plots different optimum tra- 
jectories in the (w,y) space for E = 10. Trajectories for higher 
final altitude (lower final speed) are started with lower X^. The vari- 
ation of X has been presented in Figure 4. Let us assume that the 
upper boxind of X is X^^^^ = 2. 7 5. Then from Figure 4, all trajecto- 
ries with X^ < 2.75 are pure variable X trajectories since the condition 
X = X„_^ is never reached. On the other hand, to generate the 
remaining optimum trajectories, we must start with X = for a 
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w 



Figure 7. Optimum Trajectories for Pull-Up Maneuver 


certain distance and then switch to variable X. The integration 
starts with the state equations only using X = then at. a certain 

point called the switching point with the state (w^, Ug , "Vg), we use 
the variational equations , that is the state equations and the equations 
for X and G, as before and continue the integration until X^ = 0. We 
notice that in this example X > X, , the initial derivative of X 
(dX/dy)^ with X^ = is negative. To generate the family of opti- 

mum trajectories, we can switch at any point where dX/dy is negative. 
But to solve a particular problem with a prescribed W£ or u^, the 
switching point has to be foTind such that the final condition w = w^ or 
u = Uf is satisfied. 
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Next, we consider the same case of Figure 4 with E'*' = 10, 

but now we have X = 2.0. Then all traiectories must start with 

max ■' 

X = X . Since X < X, in this example, dX/dy is positive 
max max b t- • / r f 

initially. The constant X = X subarc must continue for a certain 
■' max 

distance until dX/dy < 0, which occurs after passage through the 
lowest point in this example. Thus all the switches occur along the 
ascending arc, with the constant X subarc longer for higher altitude 
gain (smaller final speed) trajectory. 

To give an explicit example, we solve the problem for the 
following initial and end conditions , 

w. = . 5 , u* = o 5 , V. = - 1/E '' with E* = 10; 
111 

(6.16) 

V7V. = .7 , w. = maximum , = free 

f 1 f f 

The physical trajectories are plotted in Figure 8. For the trajectory 

without lift constraint, it is found that X^ = 2. 628314 leading to a 

final value w^ = 0. 53457 corresponding to a gain in altitude of h^ - h^ 

= 478.81 meters. If the constraint X = 2.55 > X. is enforced, 

max b 

the trajectory starts with X = X until w = .49935, and switches 
to variable X. The switch occurs during the descending phase. The 
final altitude is w^ = 0. 53450 and corresponding to a gain in altitude 
of h^ - h^ = 477.87 meters. On the other hand, with the constraint 
X = 2.0 <X, , the switch occurs at w = .50877, at a point along 
the ascending arc. We obtain W£ = .53402 which corresponds to a 
gain in altitude of h^ - h. 


= 471.44 meters. The variation of the 
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igure 8. Pull-Up Maneuver with Bounded Trajectories in the 
Physical Space p (h - h^^) vs y, and Variation of the 
Optimum Bounded Lift Control 
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normalized lift coefficient X for these three trajectories is also 
presented in Figure 8. 


6o 3 Maximum Range 

In this case, it is proposed to find the optimum lift control X 
to glide the vehicle from the initial condition until the final condition 


w = , u = u^ , V = 


(6. 17) 


such that the final range y£ is maximized. Since is not free, 

Cj ^ 0, and the Hamiltonian integral Eq. (5. 7) is inoperative in our 
formulation. We still have the same differential system, that is, 
Eqs. (5,4) and (5, 9), the difference here is that we have two arbi- 
trary parameters X^ and G^. The differential system is integrated 
with a set of guessed values X^ and until the prescribed final 
altitude w = w^ is reached. The other two prescribed final values u£ 
and Vj are used to adjust the values of X^ and G^. If the final angle is 
free, the condition on Y£ is replaced by the transversality condition 


X^ = 0. 


The advantage of using the variables X and G to replace the 
adjoint variables is that their numerical values are nearly constant. 
This is because in glide for maximum range, both Y and u vary 
slowly so that dY/dy 0 and du/dy 0, and we have 

. w cos Y 

A. — ~ 

^ (6.18) 


u(l + X^) 


2 E* w cos Y 


-tanY = 



38 


Hence 

- 

In this so-called steady state approximation, the range is maximized 
by selecting the lift coefficient X to minimize the glide angle - Y. 

This leads to the selection X = 1, that is, to glide at maximum lift- 
to-drag ratio. Then we have the minimum glide angle 

- tanY = (6.20) 

E* 

Of course, this solution is only approximate. The real optimum 
solution is obtained by a lift modulation. Nevertheless, the steady 
state solution provides an educated guess for the behavior of X and G. 
First, using X = 1 and Eq. (6.20) in the first equation of Eqs. (5.9), 
we deduce that dG/dy ai 0. This means that G is nearly constant 
during the glide. Furthermore, using X = 1 and Eq. (6. 20) in the 
second equation of Eqs. (5.9) and noticing that dX/dy 0, we have 
G 4 since u u. = . 5. In summary, the range of values for X^ is 
close to 1 and the range of values for G^ is close to 4. In other 
words, the optimum trajectories are very sensitive to the initial 
values Xj^ and G^, especially when E is large. 

The results are presented in Figures 9, 10, and 11 for the 
case of fighter aircraft with E* = 10. Each figure presents several 
optimum trajectories with different final altitude. To restrict the 
plot to a one-parameter family of trajectories, we impose the 
condition U£ = w£ at the final instant. For each trajectory, that is. 
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for each prescribed final altitude W£, the corresponding initial value 
of the nornaalized lift coefficient is also labelled in the figures. 

For comparison, in each figure we plot in a dashed line the steady 
state trajectory, that is, the trajectory generated by using X = 1. 

Figure 9 gives the variation of X as a function of y for 
different altitude drops. It is seen that, for large altitude drop, 
optimum glide is effected at near maximum lift-to-drag ratio, i. e. , 

X 1, except for the initial phase and the final phase. Also, we 
assume that is free, thus X^ = 0. 

Figure 10 gives the variation of -V as a function of y. For 
large altitude drop, it is steadily increasing at a very slow rate 
except for the initial and the final phases. Hence, as an approxima- 
tion, along this portion of the optimum trajectory Y is nearly constant. 

Figure 11 gives the variation of the dimensionless dynamic 
pressure t| = u/w. For large altitude drop, it is nearly constant and 
slightly less than unity during the main portion of the glide. To find 
this near constant value, we take X = 1 in Eqs. (6. 18) and have 

E* 

T] = . ( 6 . 21 ) 

VI + E*^ 

For E = 10, this value is = .99504 and is slightly less than the 
optimum value of "n which is near t) = . 9965. 

Concerning the actual performance, namely the maximized 
range, the X = 1 trajectory gives a good approximation for large 
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Figure 9. Variation of the Optimum X for Glide with Maximum Range 
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Figure 10. Variation of the Flight Path Angle for Glide with Maximum Range 
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CO 



Figure 11. Variation of the Dimensionless Dynarnic Pressure for Glide with Maximum Range 
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altitude drop. Its range is within 1% of the optimum range. A better 
approximation is to use constant dynamic pressure glide. This also 
has the advantage of flying with constant indicated speed. The value 
of p is given approximately by Eq. (6.21), and is purely a function of 
the maximum lift-to-drag ratio E"^. 

In contrast with long range glide, the optimum glide for small 
altitude drop is not close to the glide with X = 1. This is shown in 
Figure 12 where again the dashed line represents the X = 1 trajectory. 
The short range problem is closely related to the problem of a pull-up 
maneuver with prescribed range. In this respect, we have the final 



Figure 12. Variation of the Optimum Lift Coefficient for 
Short Range (Small Altitude Drop) Glide 
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condition either 


Wf , ~ prescribed , = maximum (6, 22) 


or 


Uf , , y£ = prescribed , = maximum (6.23) 


We terminate the integration at y = y^, and use the other two pre- 
scribed final values to find the two initial values X. and G. . If Y^. is 

Ilf 

free, the final condition in Y is replaced by the condition X^ = 0. 


6. 4 Linearized Singular Lift Control for Maximum Range 

An inspection of the data presented in Figure 9 clearly shows 
the difficulty encountered in the numerical computation. More 
explicitly, for large altitude drop, the initial value X^ has to be found 
with great accuracy for the final condition to be identically satisfied. 

It is seen that, except for the initial maneuver and the final maneuver, 
the lift control nearly follows the same line. This line can be con- 
sidered as a singular arc familiar to the problem in which the control 
is linear. To reduce the computation work, if this singular arc can 
be fo\ind, one can follow the line until near the end and then compute 
separately the last arc where again, the control undergoes drastic 
change. 

In general, let us consider an optimum control problem with 
the Hamiltonian 


H = H (p, X, u, t) 


(6.24) 
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In general, the solution of Eq. (6.26) provides the optimum control 

u* = u*(p, X, t) (6.27) 

This control is of the Euler- Lagrange type and u* can be expressed 
explicitly in terms of the state vector x and possibly the time t and 
some constants of integration provided that the adjoint vector p can be 
expressed in terms of the same variables. This, in turn, requires 
the analytical integration of the equations for the adjoint vector p. 
But,unfortunately , for most realistic cases it is not possible. 

Now, let us assume that we know an approximate law for the 
optimum control, say 

u"^ cii Uq (x, t) (6. 28) 

Then by Taylor’s series, we can expand the maximized Hamiltonian 
near the value u = Uq to have 

H* = H (p, X, Uo, t) + (|“)^ (u* - uo) + • • ' (6, 29) 

If Uo is near the optimum value, the difference e = u - Ug is small 
and, by retaining only the first order we have the approximate H 
which is now linear in u. Again, for this linearized problem, the 
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optimum solution is either u = u . oru = u or a variable u when 
^ min max 

the coefficient of the linear control, called the switching function, 
vanishes identically. That is 

Here, since we have assumed that u is not on the boundary, it is of 
the variable type. We have Eq. (6.30) which provides a relation 
between the state variables and the adjoint variables. This relation 
is exact in the linearized problem, but is approximate in the original 
non-linear problem. The accuracy of this approximation is of the 
order of £ . In the linearized problem, the Eq. (6. 30) is valid as long 
as the control is of the interior type. Hence, we can take its deriva- 
tive with respect to the independent variable t, to generate another 
relation between x and p. It is known that we can take the derivative 
successively until. the linear control appears for the first time, with 
an even derivative. The linear control can then be deduced explicitly. 
Then in the case where it can be expressed explicitly in terms of the 
state variables by using the additional relations obtained, we have an 
approximate but explicit law for the optimum control. 

As an example for our present case here, from the steady 
state and the numerical analyses we have found a good approximation 
for the lift control \. It is Xq = 1. We shall call this the zeroth order 
solution. By applying the linearizing technique on this problem, we 
can obtain the approximate law for the control 
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X, = 


w cos V 
u 


(6,31) 


We call this the first order solution. Then if we apply the linearizing 
technique once again, based on the first order solution which is a 
better approximation than the zeroth order, we finally have 

wcosV A 


Xe = 


u 


B 


(6.32) 


where 

A = E*^ u^ [ 6tan^ Y - 4u (1 - tan^V) + Zu^tan^V + (1+u) C] 

- 2E*u(3+u) (2 + C) wsinY + 2 C (2 + C) w^cos^Y 

B = E*2 u^ [2 (1 +u)tan^Y - 4u - C] - 2 E*u (4 + C) w sin Y 

+ 8 w^ cos^ Y 


(6.33) 


with 


C = 1 - 


u 


w^ cos^ Y 


(6,34) 


Eq. (6. 32) gives the explicit second order solution for the lift control. 
The details of the derivation of Xx and Xj will be given in Appendix B. 
For the first order solution, from the third equation of (5.4) we see 
that using the near optimum law (6. 31), the flight path angle is 
maintained constant, a fact which can be observed in Figure 10. It 
is an improved approximation as long as the zeroth order solution 
Xq = 1 is accurate. Then for the second order solution (6.32), we 
have tested it numerically, and it gives excellent results. Using the 
initial values of the No. 2 trajectory in Figures 9, 10, and 11, we 
start the integration optimally. Then at y = 2. 4 it is switched to the 
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explicit control law (6.32). The integration keeps on going, and the 
trajectory is generated by using this approximate control law. As 
compared with the nearly linear portion of the optimum trajectory 
No. 4 in Figures 9, 10, and 11, the two lift coefficients, approximate 
and optimum, agree to four significant digits, and the two trajec- 
tories generated are identical. 



CHAPTER 7 


OPTIMAL TRAJECTORIES FOR SPHERICAL EARTH 

As in the case of the flat earth, we shall consider two types 
of optimum maneuvers. The first type is the pull-up maneuver, 
and the second one is of the gliding type. The optimum trajectory 
can be initiated from the top of the atmosphere. In some cases the 
pull-up maneuver gives the skip trajectory. Since the state equations 
we have derived for the spherical planet case, the Eqs. (4. 2), are 
exact, they are also valid for' the Keplerian motion of the vehicle 
after skipping out of the planetary atmosphere. In the other cases, 
the vehicle may reenter the planetary atmosphere after a coasting 
flight to initiate a new skip trajectory until effective entry at low 
speed. We shall consider both cases. The computation is done with 
the value k^ = 900 for the earth atmosphere. For the maximum lift- 
to-drag ratio, a reasonable value E '' = 3 is considered since the flight 
is effected at high speed. Again, the maximum lift-to-drag ratio 
trajectory, that is, the \ = 1 trajectory, is used for comparison in 
the gliding type optimum trajectory. 
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7, 1 Pull-Up Maneuver at Moderate Altitude 

This is the same problem as discussed in Section 6, 1. The 
differential system consists of the Eqs, (4.2) and the Eqs. (4.7). 
The initial condition is 


0 = 0, Z = Z. , v=v^, y = 


(7,1) 


It is proposed to find the optimum lift control to bring the vehicle 
from this initial condition to the final instant 0^ such that either 


Z = Z^, 


V = = maximum 


(7,2) 


or 


V = , Z = Z^ = minimum 


(7,3) 


The final range 0^ is assumed to be free and hence Ci = 0, Since 
the second equation of Eqs, (4.7) can be deleted, the only arbitrary 
parameter is and F, can be obtained from Eq. (4, 9). If the final 
flight path angle y^ is prescribed, it is used to find the required 
initial value \^. If is free, we have the transver sality condition 

Xf = 0, 

The problem considered here involves relatively low speed 
and altitude, and we shall take the initial values as 


Z.=.S, v. = .15, v. = -^ 


(7.4) 


At high altitude where Z. 0, and v. is of the order of orbital speed. 


V. 1 , this type of maneuver leads the vehicle to skip out of the 
atmosphere. This case will be analyzed in detail later. 
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Again, the different solid lines are the optimum trajectories leading 
to the final boundary plotted in dashed line. 

The variation of the normalized lift coefficient X as a function 
of the flight path angle y for different trajectories is presented in 
Figure 14. The behavior is the same as in the flat earth case, but 
the values of X at the lowest point, Y = 0, are not so nearly the same. 
The difference is more due to the fact that the value of the maximum 
lift-to-drag ratio E used is relatively low rather than due to the 
additional centrifugal acceleration term which is included in the 
spherical planet equations. 

7, 2 Keplerian Motion Following a Skip Maneuver 

In a skip trajectory, the vehicle enters the atmosphere at 
very high altitude with a speed at orbital magnitude and uses its 



Figure 14. Variation of the Optimum Lift Coefficient for Pull-Up 
Maneuver at Low Speed over a Spherical Earth 
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lifting capability to negotiate a turn. It is then ejected from the 
atmosphere. This maneuver is depicted in Figure 15. 

The skip maneuver is an important maneuver. It can be used to 
achieve maximum range or to assist a climb to orbital altitude with 
maximum residual speed, hence minimizing the required character- 
istic velocity for orbit insertion. In the three-dimensional maneuver, 



Figure 15. Geometry of a Skip Trajectory 
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it can be used to change the orbital plane. The first order solution 
for constant lift-to-drag ratio with the centrifugal and gravity accel- 
erations neglected has been obtained in the classical literature [l] . 
An accurate second order solution for constant lift-to-drag ratio 
skip trajectory has also been obtained [?] , In this chapter, we shall 
analyze the optimum solution with lift modulation using the exact 
equations. 

The equations we have derived are valid for flight in the 
vacuum by taking the limit Z 0. But to initiate atmospheric flight 
we must start with some nonzero initial value We shall adopt 

the convention that atmospheric entry is initiated when the accelera- 
tion due to atmospheric lift is equal to a certain small fraction of the 
gravity acceleration. From the definitions of Z and v in (2.7), the 
dimensionless acceleration due to a lift force with C, = C *, is 

J-I Li 



For the earth atmosphere, pr = 900. Taking a/g = .015, i, e, , 1.5%, 
with an initial speed equal to the orbital speed, = 1 , we have 
Z^ = . 0005. We shall use this value as the value of Z at the top of 
the sensible atmosphere. For higher altitude with Z < Z^, the flight 
is considered as in the vacuum and Kepler ian motion applies. 

As shown in Figure 15, the initial point (r^, V^, V^) is con- 
sidered as the entry point, and the final point (r^ = r^, V^, Y^) is 
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considered as the exit point. Between the two points is the 
atmospheric skip trajectory, while beyond the exit point the flight 
is in the vacuum. Once in the vacuum, the vehicle climbs to the 
highest point (r , V , V = 0), the apogee of the Keplerian orbit. 
Because of the obvious symmetry, the range angle ^ between the 
exit point and the apogee is half of the range angle for the coasting 
portion of the trajectory in vacuum. We shall be concerned with 
the maximizing of either the apogee distance r , or the apogee speed 
Vg^, or the coasting range angle 2^. Hence, it is necessary to express 
these elements in terms of the variables at the exit point where 
atmospheric flight terminates. These relations can easily be 
obtained by using the classical Keplerian equations. However, we 
shall derive the pertinent equations from the general equations (4. 2). 

With Z “*■ 0, and using the equation for the variation of the 
radial distance to replace the first equation of (4.2) since it is 
inoperative, we have 



^=-(2-v)tanV 


d~Y _ (1 - v) 

d9 V 


From the second and third equations of (7. 7), 

dv V (2 - v) 


dY 


1 - V 


tan Y 


(7.7) 


(7.8) 
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Upon integrating this equation, we have 

V (2 - v) cos^Y = (1 - e^) 


(7.9) 


where the right hand side represents a constant of integration. 
Next, from the first and second equations of (7.7), 


dr _ -r 
dv 2 - V 


(7,10) 


Its integration gives 

f = 2-v (7,11) 

cL 

where a is another constant of integration. Returning to the defini- 
tion of V, V = r V^/p , it is easily seen that Eq. (7, 11) expresses the 
conservation of energy and a is the semimajor axis of the Keplerian 
orbit. Furthermore, combining the two integrals (7.9) and (7. 11) 
and again using the definition of v, we have 


r^V^cos^V = pp (7.12) 

where 

p = a(l-e^) (7.13) 

Equation (7, 12) expresses the conservation of angular momentum, 
and it is now clear that e is the eccentricity of the orbit while p is 
the semilatus rectum. Now, consider the derivative 



p dr 
r^ d0 


^ tan Y 
r 


(7.14) 


where we have used the first equation of (7. 7). By taking the 
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derivative again and simplifying, we have 

^ (?) = (?)[- (7.15, 

Using the definition of v and the integral (7. 12) in Eq, (7, 15), we 
finally have the differential equation for the orbit 

The general solution of this equation is 

^ = 1 + Acos e + B sin 0 (7.17) 

r 

where A and B are constants of integration. Starting the angular 
variable at the perigee, 0 = 0, dr/d0 = 0, r = a (1 - e), and we obtain 
the polar equation of the orbit 

^ = 1 + ecos0 (7.18) 

Hence, we have derived the classical equations for Keplerian motion 
from our general equations (4.2). With these equations, we can 
deduce the performance indices for optimization in the following 
sections, 

7, 3 Skip Trajectory for Maximum Final Speed 

Again, we assume that 0^ is free. Thus Cj = 0, and we can 
use Eq. (4, 9), It is a one-parameter problem. Referring to Figure 
15, the vehicle enters the atmosphere at the initial point with the 
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initial condition 

(Zi , , V.) = (. 0005 , 1,0 , variable) (7. 19) 

It is proposed to find the optimum lift modulation such that at the 
exit point 

= Zj , ~ » ■'^£ = maximum (7, 20) 

Since is free, we again have = 0. We integrate the Eqs, (4, 2) 
and the first and third equations of (4. 7), from the initial state, with 
a guessed X^ and a solved from Eq. (4. 9), Then we use the con- 
dition X^ = 0 to find the correct value of X. , 
f 1 

The variation of the optimum lift coefficient as a function of 
the speed ratio V/V^ is presented in Figure 16 for several initial 
flight path angles It is clear that less negative gives higher 
final speed. For all the trajectories computed, the optimum lift 
coefficient slightly increases at the beginning and then decreases 
continuously to the final value X^ = 0. 

7. 4 Skip Trajectory for Maximum Apogee Altitude 

In this problem, it is proposed to use optimum lift modulation 
to bring the vehicle to the exit point such that subsequent climb 
in the vacuum leads to a maximum height. Since 0£ is free, Cx = 0, 
and therefore the only parameter is X^ . From Eq. (7.11) we have 

r 2 - V 

— =T ^ (7.21) 

■■f ^ - ''f 






60 


Since J is a function of and we have at the final time 


3J 

3Vf 


3J 

9Y. 


(7.24) 


Upon using the relation of Eq. (6. 5), we have the following trans- 
versality condition 


X, = 


E* (2 - v^)^ sin 'Y£ cos Y£ 


^ 2 [l - (2 - V£) cos^ Yf + -^1 - (2 - Vf) V£ cos^ Y£ ] 


(7.25) 


This condition is used to find the initial value X^ for the optimum 
trajectory. Finally, the corresponding max (r /r,) can be obtained 

3> X 

from Eq. (7. 23), 

This problem has been solved, and we have the following 

r e suits 


= .0005 , 

V. = 1 . 0 , 

o 

CO 

1 

II 


1 

1 

= . 0005 , 

V£ = . 377 , 

Y£ = 43.36° 
/ r \ 

- .70225 , 

= 2.0440 6 , 

max ( — ) = 1. 12308 


f 

ViJ 


In this flight program, the initial lift coefficient is negative. It 
appears that the optimum trajectory starts with a plunge toward the 
dense atmosphere with a slight increase in the speed, and then uses 
the lift to rotate the velocity vector upward with a relatively high exit 
angle, to achieve the absolute maximum apogee height. This 
maneuver is purely an academic exercise. It incurs excessively 
high acceleration. The value of Z at the bottom of the flight path is 
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= 4 . 7795 , This represents a dip into the atmosphere with a 
distance Ah = (1/p) log (Z^/Z^) = 65,641. 5 meters. In practice, we 
shall have the following problem. 


7. 5 Skip Trajectory for Maximum Apogee Altitude with Prescribed 
Apogee Speed 

This is a more realistic formulation of the previous problem. 
Specifically, we seek to maximize the apogee altitude while prescrib- 
ing a residual apogee speed V . 

3 > 

L.et 



V ^ 
— a 

V = — ; — - 

a V-lr^ 


(7.26) 


Hence, we minimize J = 1/R with a prescribed v . From Eqs. 
(7.21) and (7.26), we have 


1 1 I 1 — 1 

R 2 a 2 £ 


(7.27) 


From this relation it is obvious that to minimize 1/R we simply 

maximize the final speed v^. But this time, besides maximizing v^ , 

the prescribed V (or v ) must also be achieved. Since v = r V^/u 
^ a a a a a 

= R V , by using this relation and Eq. (7.27) for 1/RinEq. (7,22), 

cL 

we have 


V 


T-~ - (1 + T V - \ V )^ = 0 

V ros^Yr ^ a ^ f' 


(7.28) 


This is the final condition to be satisfied so that the prescribed v_ 


can be achieved. The procedure to solve this problem is the same 
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as that in the preceding section, but here we use Eq, (7. 28) to 
search for the correct initial value 

1 

As this is a one-parameter problem, we can obtain the 
totality of solutions simply by varying the parameter Then at the 

end of the integration where = Z^, Eq. (7.28) is used to solve for 
V and Eq. (7.27) for R. Figure 17 presents the solution for several 

3 > 

values of Y. . For each value of Y. there is an absolute maximum 
1 1 

apogee distance corresponding to the problem solved in the preceding 
section. For any other prescribed v which is different from this 
point, the maximized apogee distance is lower. The case of v =0 

cL 

corresponds to vertical ascent in a vacuum, and hence for a tra- 
jectory leading to Y^ = tt/2. 

7. 6 Skip Trajectory for Max:imum Apogee Speed with Prescribed 
Apogee Altitude 

This is a trajectory with practical importance. It is 

proposed to use optimum lift modulation to bring the vehicle to the 

exit point such that the subsequent ascent in the vacuum will lead 

the vehicle to a prescribed apogee altitude r with a maximized 

sl 

residual speed V . Clearly, this leads to minimizing the character- 

Si 

istic velocity AV for orbit insertion. 

By eliminating v^ between Eqs. (7. 27) and (7.28) we have 

(2 - v^) - 2 R + v^ cos^ Y^ = 0 


(7.29) 
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For a prescribed r , the ratio R is specified. The procedure to 

3l 

obtain the optimum solution is the same as in Section 7,4, except 

that in here Eq, (7. 29) is used to adjust the initial value X. , For 

the totality of the solutions, it is exactly the same as has been plotted 

in Figure 17, But this time, the value r /r. is prescribed while the 

a I ■ 

corresponding value v is maximized. We notice that there exists a 
range of that gives two trajectories both satisfying the neces- 

sary condition for optimality. The optimum trajectory is the one 
corresponding to higher value of v . 

7, 7 Skip Trajectory for Maximum Coasting Range 

Again, we refer to Figure 15, For the initial condition we 
are still using Eq, (7, 19), In this problem, it is proposed to find 
the optimum lift control to negotiate a skip trajectory such that after 
its exit from the atmosphere, the vehicle coasts ballistically in the 
vacuum to achieve a maximum coasting range 2^. We first solve 
this problem by assuming that the final value 0^ at the exit point is 
free, hence Cj = 0. This is suggested by the fact that at orbital 
speed with small value of Y., the coasting range 2^ is significantly 
larger and more sensitive to change than the atmospheric skip range 
(0^ - 0^), The next case to be addressed is the maximization of the 
total range from the initial point, (0^ - 0.) + 2^ , 

From Figure 15, it is seen that ^ = tr - 0£, Therefore, we 


obtain from Eq. (7. 18) 




65 



(7.30) 


By writing the Eq. (7„ 9) at the exit point and solving for e, we have 


e 


- (2 - v^) v^cos^Y^ 


(7,31) 


Then by using Eq, (7. 31) and the Eq. (7. 11) at the exit point in the 
relation (7, 13), it gives 

p = r^v^cos^Y^ (7.32) 

Upon substituting Eqs, (7.31) and (7.32) into Eq. (7.30), it becomes 

1 - V£Cos^Yf 

cos i = — -■ : (7. 33) 

yjl - (2 - v^) v^cos^ Y^ 


For the first case we maximize 2^, It is equivalent to minimizing 
cos ^ and thus J = cos Since J is a function of the two final 
variables v^ and Y^ , we again have the relations (7, 24). The trans- 
versality condition is then 


E* [l 




- v^ - tan^Y^] 
2 tan Y r 


(7,34) 


This is the condition used to search for the exact value The 

1 

initial state used and the results obtained are 


Z. = .0005 , V. = 1.0 , Y. = - 4° 

1 1 .1 

Z^=. 0005 , v^=. 87475 , Y^ = 6. 02° 
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= .2925 

1 


0 e, = . 17 646 
f 1 




max 


1, 18958 


To show the optimality character of this trajectory, we integrate the 
state equations (4, 2) using a constant lift coefficient, \ = constant. 
The best constant X which gives the maximum coasting range is fo-und 
to be X = 1.024, and the results are 


Z^=.0005 , v^=. 90876 , \^ = 3.58° 


- 0. = . 20633 
f 1 




It shows that using optimum lift modulation we have an improvement 
of 10.41% in the coasting range as compared to the best solution 
obtained with a constant lift coefficient. 

We now solve the second case, in which we maximize the 
total range from the initial point to the end of the coasting flight. 
That is, we maximize the following performance index 

[ 1 - Vr cos^ Vf 

— • 

\ 1 - (2 - v^) v^cos^Y^ 


(7,35) 


This time, the final range is not free and hence Cj ^ 0, There are 
two parameters to be found, X^ and F^, Actually Cj is equal to Pq , 
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and from. Eq, (7. 35) we have, for a maximization problem, 




(7.36) 


Thus Cx = 1, Furthermore, we also have 


8J 
~ 3v, 


aj 

av. 


(7.37) 


Upon using the relations { 1 . 31 ) in Eq. (6.5), it gives exactly the 
same as given by Eq. (7.34). Therefore, Eq. (7.34) is also a 
trans ver sality condition for this case. We need one more transvers- 
ality condition because this case has two parameters. It comes from 
the Hamiltonian integral (4. 8) at the final time. With Cx = 1 and the 
given in Eq, (7.37), we finally have 


f 


kZ v (l-X=) (l-v)\ V 

+ g? + (l - Y + VjFjjtanVj = 0 (7.38) 


E'^ cos V , 


The problem is solved and this time it is found that 


= ,0005 , v^ = . 88101, V^ = 5.63° 

= .57921 , F. = 3. 6873 
1 1 

9^ - e. = . 18173 
f 1 

2 § = 1 . 18692 

The total range obtained is J = (0^ - 0^) + 2^ = 1.36865, which is 
slightly higher than the total range J = 1. 36604 of the first case 
where only the coasting range 2§ is maximized. 
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For all the skip trajectories solved above, it has been 
assumed that beyond the altitude Z = . 0005, the flight is in a 
vacuum. The corresponding altitude is given, through the atmos- 
pheric density, by 


2m 



= .0005 


(7.39) 


For most vehicle characteristics, the resulting altitude is generally 
high enough such that beyond this altitude the subsequent trajectory 
is practically Keplerian. For better accuracy, one can take a 
smaller Z. The computational procedure remains unchanged. 


7. 8 Glide with Maximum Range 

The maximum range obtained previously concerns the range 
with one skip. We now generate the optimum control to maximize 
the total range for a descent from an initial altitude Z^ to a final 
altitude Z^ . The problem is first solved for the case of a relatively 
low initial altitude, A reasonable set of initial values is 

(Z. , V. , Y.) = (.5, .15, - (7.40) 

1 1 

with again E '' = 3. This can be considered as the gliding flight 
following a ballistic entry of a shuttle vehicle. The vehicle enters 
the earth atmosphere at the reentry altitude Z^ 0, with a speed 
v^ ^ 1 and a certain reentry angle Y^. Then at the end of the 
ballistic phase, the vehicle rotates to reduce the angle of attack. 
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hence generating a lifting force and using lift modulation to glide to 
a final altitude with a mciximum range. The case of gliding from 
the entry point will be analyzed in the last part of this section. 

The numerical computation is carried out exactly as in the 
flat earth case. We integrate the state equations, Eqs. (4.2), the 
first equation of (4.7) for F, and the third equation of (4.7) for \, 
from the initial state (7.40) and two guessed values and , At a 
prescribed final altitude the integration is terminated, and the other 
two prescribed final values v^ and are used for adjusting the F^ 
and Xj . If Y£ is not prescribed, then the condition on V£ is replaced 
by the condition X£ = 0. In order to generate a one-parameter family 
of optimum trajectories, we impose the final condition 

V?rZ^V£=l (7.41) 

Physically, this means that the final acceleration due to a lift force 
with is equal to the gravity acceleration. For each pre- 

scribed Z£ , the corresponding v^ is obtained from this condition. 

Figure 18 presents the variation of the optimum lift 
coefficient. It is seen that X oscillates about the value of unity and 
tends to this value near the end of a long range glide which corre- 
sponds to a large altitude drop. Figure 19 presents the variation of 
the flight path angle while the variations of the altitude and the speed 
are depicted in Figure 20. 






P(h-h, 
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.05 ,10 .15 .20 .25 

Figure 20. Variation of the Altitude and the Speed for Low Altitude 
Maximum Range Glide over a Spherical Earth 
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For the case of glide starting from the entry point, a typical 
initial condition is 

(Zi , Vi . Yi) = (. 0005 . 1.0 , -4°) (7,42) 

The final condition to be satisfied is 

Z = , V = , Y^ = free (7, 43) 

The variation of the altitude as a function of the range for the optirtial 
trajectory is plotted in Figure 21 as a solid line. The variation of 
the flight path angle is plotted in Figure 22. Finally, Figure 23 
presents the variations of the optimum lift control and the speed. 
Again, the optimum lift control oscillates and tends to the lift control 
for maximum lift-to-drag ratio, X = 1. 

In both cases above, the trajectory generated by using 
maximum lift-to-drag ratio, X = 1, is plotted in the dashed line for 
comparison. Besides an improvement in the range of about 2%, the 
oscillation in altitude along the optimum trajectory is less severe. 

We can also see a more desirable behavior of the flight path angle 
along the optimum trajectory. It also yields a more smooth variation 


in the deceleration. 



,0005 
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igure Zl, Variation of the Altitude for Maximum Range Glide from the Reentry Point 
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Figure Z3„ Variation of the Optimum X and the Speed for Maximum Range 
Glide from the Reentry Point 



CHAPTER 8 


THREE-DIMENSIONAL FLIGHT 

Two elements are of interest in three-dimensional flight; one 
is the maximum cross range, the other is the footprint. In this 
chapter, we shall discuss the procedure to solve the absolute maxi- 
mum cross range. We will see that it is a three-parameter problem. 
Then by using the equilibrium glide condition as a simplifying device, 
we shall compute the footprint of a gliding entry vehicle on the surface 
of a planet. A technique of coordinate rotation is used to make 
the iteration much more effective. 

8. 1 Maximum Cross Range 

It is proposed to find the lift and bank modulation to maximize 

the final latitude <()^ while the final longitudinal range 0^ is free. For 

an initially circular orbit, if the position of departure is free, the 

reachable domain will then be a zone between the latitudes -i 

^mcix 

and +cb . If d) = tt/Z, the reachable domain is the entire 
mcix max 

surface of the planet. 

Since the final arc length is free, we have Cq = 0 in Eq. 
(3.6). The final condition in the state variables will be 
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Z = , V = , y = y^ = free , 0 = 0^ = free , 


4' = = free , 4> = = maximum 


( 8 . 1 ) 


The Eqs. (3.7) at the final time then can be written as 


Ps ^ c, = 0 


p , = - cos 4>^ (C2 sin 0^ - C3 cos 0^) = 0 


( 8 . 2 ) 


= C 2 C 0 S 8^ + C3sin0^ = 1 


Solving for the constants of integration C2 and C3 we obtain the 
solutions for p^ and p^ , 


p , = cos <f> sin (0 . - 9) 

4 * 1 


p^ = cos (0^ - 9) 


We also have, since 7^ is free. 


(8.3) 


Qf = = 0 


(8.4) 


The Hamiltonian integral (3. 11) becomes 

kZ 


- P tan 7 - N 




(2 - v) 

* + tan 7 

E*cos7 ^ 


(1 - v)Q 

V 


4 N cos 7 


E^kZ ly 

+ ~m ( + — 2 — / + p , sin Jj - p , cos 4 tan <b = 0 (8.5) 

‘ cos‘7/ 4 '^4 ^ ^ 


Thus for the specified final condition (8. 1), the procedure to obtain 
the optimum solution is as follows. Starting from a certain initial 


state, say 
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(Zi, V., V., 0., 4>., 4^.) = (.0005, 1.0, - 40 , 0, 0, 0) (8.6) 


with a given value of we integrate the state equations (2.8) along 

with the adjoint equations (3. 9), using the control law (3. 10) and 

the Eqs. (8.3) for and p^. There are three parameters, namely, 

the final longitudinal range 0^ and two of the three initial values P. , 

N. , and Q. since one of them can be obtained from the Hamiltonian 
1 1 

integral (8,5). These parameters are to be selected such that when 
the integration is stopped at 0 = 0^, the two prescribed final values 
and v^ and one trans versality condition (8. 4) are all satisfied. 

The resulting trajectory will be the optimum trajectory for maximum 
cross range, 

A simplification can be made by using the so-called equilib- 
rium glide condition, assuming that the glide angle is small and 
stays nearly constant. This is expressed as 


y ^ 0 



(8.7) 


By substituting into the equation for y in (2.8), we have 


kZ = 


1 - V 

X V cos O' 


( 8 . 8 ) 


This equation is used to evaluate the altitude Z. Thus we have the 
following reduced set of state equations 
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dv _ (1 + X^) (1 - v) 

ds E*X cos cr 


d9 _ cos 4^ 
ds cos 4> 


(8.9) 


d4> . I 

~r^ = sin ip 
ds 


ds 


(1 - v) 

V 


tan O’ - cos 4* tan 4> 


The Hamiltonian of the reduced problem is 


H = - 


(1 + X^) (1 - v) cos 4 

P + p„ r + p, simp 

E^X cos cr ^9 cos tp 


P I. tan cr - cos 4 tan cpj 




( 8 . 10 ) 


Then, it is clear that the optimum lift control is 


X = + 1 


( 8 . 11 ) 


that is the glide is effected at maximum lift-to-drag ratio. For the 
bank control, we either have 


cr = O' 


max 


or an interior bank control such that 


sin O' = 


E*p, 

2 vp^ 


4 


( 8 . 12 ) 


We notice that the integrals (3.7) are still valid for this case. Hence, 
with Cl = 0, Cq = 0 , we can write the Hamiltonian integral of (8.10) in 
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the form 


2py (1 - v) 
E* cos (T 


+ p , sin 4^ + p I 


I 


(1 - v) tan O' 


cos 4* tan d) = 0 


) = 


(8, 13) 


Using the optimal law (8. 12) to eliminate p^ and the Eqs. (8. 3) for 


p^ and p^, we have the explicit law for the bank angle 


tan cr = 


(1 - v) 


cos <j> sin (0^ - 9) 


V cos (0^ - 9) sin 4* - cos sin (ji sin (9^-9) (8,14) 


The problem is thus reduced to a one-parameter problem in the 
parameter 0^. In this formulation the stopping condition is no longer 
but the final speed , 

For numerical computation, we use the control law (8. 14) to 
integrate the full set of exact state equations (2, 8) with a guessed 
value for the final longitudinal range 9^. This value is to be adjusted 
such that, at the final time when 9 = 9£, the prescribed final condition 
V = v^ = . 001 is satisfied. The initial state used is (8. 6) except that 
the initial speed is 0.99 instead of 1.0. The purpose of this change 
is to give a defined cr value at the initial instant. The maximum 
value of the bank angle is selected to be 8 5°. Figure 24 presents 
the maximum cross range solved by using the reduced control law 
(8.14), as a function of the maximum lift-to-drag ratio E*. The 
dashed lines represent the results of the gliding trajectory with \ = 1 
and cr = 45° where the bank angle is switched to 0° when the heading 
angle 4^ reaches the limiting value 90°. The improvement in the 
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1 , cr from Eq. (8. 14) 


1 . tr = 45° and 0° 


Figure Z4, The Longitudinal and Cross Ranges as Functions of E''" 


cross range is easily seen. We can also see that for vehicles with a 
maximum lift-to-drag ratio E* greater than the value 3. 5, the maxi- 


mum cross range is larger than 90° and the reachable domain of the 


vehicle is the whole surface of the earth if it has an initially circular 
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orbit and the point of departure is not specified. For all the trajec- 
tories, the final altitudes are about the same with 30, which 

corresponds to an altitude drop of about 50 miles from the initial 
point. The final flight path angles vary with E , with larger E giving 
flatter flight path angles. For example, for E =1.5 the final flight 
path angle is = - 23. 5° for both control laws Eq. (8. 14) and cr = 45° 
and 0°j for E"*' = 3. 5, it is Y^ = - 11. 0°. Figure 25 presents the 
variation of the altitude and the speed of the trajectory generated by 
the lift control X = 1 and the bank control (8. 14) while Figure 2 6 
presents the variation of the flight path angle and the bank angle, for 
the maximiim lift-to-drag ratio E =1,5. 

8. 2 The Footprint 

As has been mentioned before, if the reentry vehicle is 

initially in a circular orbit and the position for leaving the orbit is 

not prescribed, then the reachable domain on the surface of the earth 

will be a zone between the latitude -4> and 6 . The footprint of 

max max 

a reentry vehicle is defined as the curve limiting the reachable 
domain on the surface of the earth if the reentry point is specified. 
This problem is even more complicated since we have to find the 
maximum cross range for each prescribed final longitudinal range 
0^. As the final longitudinal range is no longer free, pg = Cj ^0. 

In Eqs, (3.7), if we divide all the equations with Ci , they become 



P(h - h. 



Figure 25. Variation of the Altitude and the Speed of the Trajectory Generated by the 
Control Law X = 1 and Eq. (8. 14) 
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= 1 


^ C, C, 

7T = cos 9 + -pr sin 0 

Cl Cl Cl 


(8. 15) 


4^ Cs „y 

— = sin <p - cos 9 I — sin 0 - — cos 9 I 

Oi 'Cl Cl / 


Similarly, Eqs. (3.9) can be rewritten in the form 


d / P 


Z 


ds VCi/ cosY 


(E*\Ci; 4 (N/Ci)L\Ci/ cos^yVCi/J/ 


ds 


(«)■ 


A. /-Q 

ds \Ci 



(8. 16) 


In terms of the new variables p,/Ci , etc. , the Hamiltonian integral 

9 

(3. 11) becomes 



P \ / N 

5 r) - (c 



kZ 


L E* cos Y 


+ 


E*k Z 


4 (N/Ci) cos Y 






sin 4^ 


where again Cq = 0 since the final arc length is always free. The 


control law is , from Eqs. (3.5), 
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X 


cos O’ = 


E-(Q/Ci) 

2(N/Ci) 


X 


sin cr = 


E*(Pj,/C,) 

2 (N/Ci) cos V 


(8, 18) 


The Eqs. (8. 15) - (8. 18) conabined with the state equations (2.8) are 
the equations for solving the exact footprint of a reentry vehicle from 
a specified departure point. There are four parameters in this 
problem; they are (C 2 /C 1 ), (Cs/Ci), and two of the three initial values 
(P/Cj)^ , (N/Ci)^ , and (Q/Cj)^ since one of them can be obtained from 
the Hamiltonian integral (8. 17). Among them, one can be used as a 
scanning parameter since we want to solve the whole footprint. 

Hence, it is a three-parameter problem. For a typical example, a 
vehicle is initially at the specified point (8. 6). To find the exact 
footprint we pick a scanning parameter and guess the other three 
parameters, and start the integration of Eqs. (2.8) and (8. 16) along 
with the using of Eqs, (8. 15) and (8. 18). The three guessed 
parameters are to be adjusted such that when the integration is 
stopped at the final time with v = v^, the prescribed Z = and 7 = 7^ 
and the transversality condition 



are all satisfied. Then by varying the scanning parameter the foot- 
print can be solved. If the final flight path angle is not prescribed, 
the condition Y = will be replaced by another transversality condi- 
tion, namely (Q/Ci)^ = 0. 
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In the preceding section we have obtained hy using the 

equilibrium glide simplification. We shall use the same device in 
this section again. The control law (8. 12) is still valid in this case. 
By using it in the Hamiltonian integral (8. 10), an explicit law for the 
bank angle is solved to be 


A 

tan 0- = - 


( 8 , 20 ) 


where 


= /" I - 

V V / Cl 


cos 4^ 

B = + sin 

cos (f> Cl 


( 8 . 21 ) 


4^ - sin 41 tan (j> 


From the second and third equations of (3. 7) since Ci ^ 0, we have 


P 4 ; 

TT" = sin 4> - cos (ki sin 0 - kj cos 9) 
Cl 


= ki cos 0 + k 2 sin0 


( 8 . 22 ) 


where 


= cT • ' Q 


(8.23) 


Hence, there are two parameters ki and kz in the reduced problem, 
and it is a one-parameter problem since either ki or kj can be a 
scanning parameter. For the transversality condition since the final 
heading angle 4<£ is still free in this case, we have p^ = 0 or from 
the first equation of (8. 22), 
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sin - cos <j)^ (ki sin 0^ - k 2 cos 9^) = 0 (8. 24) 

The procedure to solve the footprint is as follows. Using the explicit 
control law (8. 20) in the full set of exact state equations (2. 8), the 
integration is started from the initial state (8. 6) with = . 99 instead 
of 1.0. For the two parameters kj and k 2 we pick kj as the scanning 
parameter and adjust k 2 such that when the specified value v = v^ is 
reached, the transversality condition (8. 24) is also satisfied. By 
varying ki and doing the same adjustment on k 2 for each value of kx , 
the whole footprint is solved. Although in this reduced problem the 
final altitude is not specified, according to the numerical results it 
is acceptable in general. 

A technique of coordinate rotation has been introduced by 
Fave [ 9 ] for a flat planet model. Its application in the spherical 
planet model enables us to use the control law (8. 14) which corre- 
sponds to Cl = 0 for solving the footprint. We shall illustrate the 
technique in the flat earth case at first, and then use it in the 
spherical earth case. In Figure 27, let M. y z be the initial coordi- 
nate axes and be an optimal trajectory leading to the final 

point on the footprint C for a given longitudinal range y^. Let 
M. y' z' be the rotated coordinate system with the axis M. y' parallel 
to the tangent of the footprint C at the point M^, Since the footprint 
is the same if the initial condition is maintained, if we use the new 



90 



Figure 27, Rotation of the Coordinate Axes, 

Flat Planet Case 

axes M. y' z' te find the point on the footprint, we have the problem 
of mclximizing z' while y' is free. As the axis My' is not known a 
priori, we have a new parameter, namely the initial heading angle 
with respect to the new axes. Besides the new parameter i|j^' , we 
also have another parameter y^' on the new axes. These two new 
parameters correspond to the parameters kj and k 2 in Eqs. (8.22), 
but, they are geometrical quantities on the rotated axes. The 
parameter i}j/ is the initial heading angle with respect to the rotated 
axes. It will be clear later on that, if we consider the upper half of 
the footprint and translate the rotated axes to the. points on the foot- 
print, we will see that the maximum longitudinal range point 
corresponds to the value 4 '.' = 90°. Then, as the new axes are moving 
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along the footprint, the 4^^' angle is decreasing from 90° to 0° and 
then to -90°. The 4^.’ = 0° corresponds to the maiximum lateral range 
point of the footprint, and the 4*/= -90° is the point where the slope of 
the footprint fails to be continuous. On the other hand, the other new 
parameter y^' is the value of y' where the maximization of z’ occurs. 
Hence, for the two new parameters we can pick 4*/ as the scanning 
parameter. For each value of 4^^’ from +90° to -90°, the y^' is 
adjusted such that the final condition is satisfied. Then, from the 
values (y^’ , z^') and the angle 4^.', we can compute the coordinates 
(y^, z^) of the resulting point on the original axes M. y z by using the 
relations 

y- = y ' cos 4^.' + z ’ sin4J.’ 
f f 1 f 1 

(8. 24) 

z, = - y ' sin 4’.* + z ' cos 4<.' 
f f 1 f 1 

By varying the 4'.* from +90° to -90°, the footprint can be obtained 
very systematically and effectively. 

For the spherical earth model, the rotation of the coordinate 
axes must be performed on the surface of a sphere, since all the 
coordinate axes must be along the great circle. The equations for 
coordinate transformation are not apparent and their derivation is 
more elaborate. Again, in Figure 28, M. 0 4> is the original coordi- 
nate system and M. is an optimal trajectory leading to a point 
on the footprint C for a given longitudinal range 0^. There is a 
tangent of the footprint C at the point Mf. At point and parallel 



Figure 28. Rotation of the Coordinate Axes, 

Spherical Planet Case 

to this line we can draw a straight line which combines with the 
center of the sphere O to decide the great circle plane for the rotated 
axis M. 9'. The rotated axis M. cji' is then on the great circle passing 
through the point and perpendicular to the great circle of axis 
M. 9'. Hence, M. 9'4>' is the rotated coordinate system for the point 
on the footprint. Referring to this new axes system, the optimal 
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trajectory has an initial heading angle 4*/, but the point is 

the absolute maximum lateral range point. Thus, with respect to 
the new axes M. 0' cj)' we are maiximizing cj)’ with 0‘ free, and the 
explicit control law (8. 14) can be used. For convenience, Eq. (8. 14) 
is rewritten on the new axes M. 0' (js' as follows 


tan O' = 


(1 - v) 


cos cj>' sin (0£ - 0* ) 


cos (Q^ - 0') sin 4^' - cos 4^' sin <f)' sin (0^' - 0') 


(8.25) 


For each value of 4^/ from +90° to -90°, we guess 9^' and integrate 

the state equations (2.8) from the initial condition (8. 6) -with v. = .99 

1 

instead of 1.0 by using the explicit control law (8.25). Then we 

adjust 9^' such that when the integration is stopped at 0' = 0^' the 

final speed v^ = .001 is satisfied. The results (0^*, c^p obtained from 

this iteration are the values on the rotated axes M. 0' cb' . The 

1 

formulas to translate them to the values referring to the original axes 
are 

tan 4>j sin 4'^^' 

tan 9 = tan0' cos 4.’ + — — 

f f 1 cos 0 ' 

(8.26) 

sin 4^ = sin 4^' cos 4.' - sin 0^' cos sin 4.’ 


These formulas are derived in Appendix C, using the spherical 
trigonometrical relations. To construct the footprint, we start from 
the value 4^' = 90° which corresponds to the maximum longitudinal 
range point of the footprint. As 4/ is decreasing from 90° to 0°, 
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which corresponds to the global maximum cross range point, the 
portion of the footprint to the right of the global point is 

obtained. The portion to the left is constructed by ranging from 

0° to -90°. Figure 29 shows the footprint for the maximum lift-to- 
drag ratio E = 1.5. The trajectories leading to the points on the 
footprint are also depicted with the corresponding values of 4^.' and 
0^* given. This technique of rotating the coordinates is not applicable 
to the short arc to the left end of the footprint beyond = - 90°. For 
all trajectories from 4^^' = 90° to -90°, the bank angle is always 
positive, that is, to the left, or zero. But for the short arc beyond 
4^i' = - 90° we have to bank the vehicle to the right at first, and then to 
the left at a certain switching point. Figure 30 presents the bank con- 
trol as a function of the longitudinal range 9 for the trajectories. The 
maximum bank angle is = 85°. For trajectories with long longi- 

tudinal ranges, the bank angle is near zero initially. It increases to 
certain value and then decreases to zero finally. For trajectories 
with short longitudinal ranges, the bank control hits the ^ 

while and then decreases to zero finally. Again, the final altitude is 
not considered. But for all the trajectories the final altitudes are 
very close to the value = 30 which is a reasonable low altitude. 




Figure 29. Footprint for 
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Figure 30„ Bank Control for Trajectories Leading to the Footprint, 



CHAPTER 9 


CONCLUSIONS 

A general solution for optimum reentry trajectories in a 
vertical plane has been presented. The three-dimensional optimal 
trajectories leading to a maximum cross range and the footprint are 
solved in a reduced problem. Unlike previous numerical studies in 
the published literature where physical data have to be specified 
numerically, here we only have to specify the most important perfor- 
mance parameter, namely the maximum lift-to-drag ratio E . The 
numerical results obtained are valid for all vehicles having the same 
maximum lift-to-drag ratio. For the other vehicles with the values 
of E'‘' around the value we have used for computation, the behaviors 
of the optimum lift control and the trajectory variables such as 
altitude, speed, and flight path angle are essentially the same. This 
advantage of having a general study is made possible by the use of the 
modified Chapman's variables and the normalized lift coefficient. 

The planetary atmosphere is assumed to be spherical and at rest, 
with locally exponential variation in its density. It is found that the 
characteristic for any atmosphere can be specified by the average 
value of the dimensionless quantity = Pr. For the numerical 
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computation, we take = 900 for the earth's atmosphere. The 
equations retain the generality and are also valid for flight in a 
vacuum. Hence, the totality of the optimum trajectory, from entry 
to landing, can be followed continuously even if at the beginning the 
vehicle skips out of the atmosphere repeatedly before effective entry 
at lower speed. 

For the planar flight case, several optimum problems for 
flight over a spherical earth are solved and the results analyzed in 
detail, especially the skip trajectory. At low altitude and low speed, it 
is more convenient to use a flat earth model. This has been achieved 
by using a canonical transformation applied to the spherical equations 
followed by a flat earth simplification. Optimum problems for flight 
over a flat earth are solved using the simplified equations. The 
optimum gliding trajectory for maximum final range, as compared to 
the maximum lift-to-drag ratio gliding trajectory, has better range 
and smaller peak deceleration and is less oscillatory. 

In three-dimensional flight, we have two more state variables, 
namely the latitude and the heading, one more control, namely the 
bank angle, and two more adjoint equations. But at the same time, 
we have two additional integrals. Hence, the real difficulty in three- 
dimensional analysis lies notin the analytical formiilation but in the 
practical computation of a two-point bo\indary value problem con- 
taining three parameters instead of two as in the planar case, A 
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simplification is thus introduced by using the so-called equilibrium 
glide condition, assuming that the glide angle is small and stays 
nearly constant. Then, by using this simplification, the footprint of 
a reentry vehicle is calculated. A technique of coordinate system 
rotation has been used, which makes the iteration much more 
effective and geometrically meaningful. 

A distinctive feature of the present formulation is that the 
equations of motion and their variational derivations are valid 
uniformly for flight in the dense layer of the atmosphere where the 
aerodynamic force is predominant and for flight in the near vacuum 
where the Newtonian gravitational force is predominant. Hence we 
can use the same equations for the investigation of the effectiveness 
of the optimum aerodynamic control at very high altitude. It is 
expected that this tenuous aerodynamic control, coupled with a thrust 
control with small magnitude, will be sufficient as optimum controls 
for the guidance of skip trajectories. 
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APPENDIX A 


Normalized Drag Polar 

Consider a generalized drag polar of the form 

where at very high Mach number, the zero-lift drag coefficient 

, the induced drag factor K, and the exponent n are assumed to 
have their constant asymptotic values. If E = is the lift-to- 

drag ratio, then 


E 



+ K C 


n-I 

L 


(A, 2) 


Hence, E is a maximum when 


—4 + (n-,1) K = 0 


(A, 3) 


This corresponds to the lift and drag coefficients 


T(n-1 




)K 


* n 




(A. 4) 


The maximum lift-to-drag ratio E‘' is, of course, 




* 


(A, 5) 


If we define the normalized lift coefficient \ as 
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^ ~ r' ^ 


(A, 6) 


then it is clear that when X = 1, and the operating point is 

at the point of maximum lift-to-drag ratio. Using (A, 4) and (A. 6) in 
(A. 1), we have 


= 


'D (n-1) 
Considering (A, 5) we obtain 


[(n-l) + X^] 


(A. 7) 


Cj3 = Cj^'f(X) 


(A, 8) 


where 


f(X) = 


(n - 1 ) + X 


n 


n 


(A. 9) 


For the case of a parabolic drag polar, n = 2, we have 
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APPENDIX B 


Derivation of the Equations (6.31) and ( 6 . 32) 

The exact Hamiltonian for the flat planet case is 


up (1 + X.^) 

wp tanY - rr^ir^ 

w E w cos Y 


Py 

- 2 p tan Y + 

u wcosY u 


-~r (B.i) 


We use the approximate solution 


Ko = 1 


(B,2) 


to linearize Then 




\=\r 


ax 


X=Xr 


(X - Xq) + 


(B,3) 


By retaining only the first order term, we have ' 


*" * A, 

= w p tan Y - 2 p tan Y - — + 77 ^ 

^w ^u u E'wcosY 


(E^p^-ZupJ (B.4) 


From the linearized Hamiltonian (B.4), we can derive the corre- 
sponding linearized state and adjoint equations. They are 


= wtanY 


dw 
dy 

du _ 

dy E* w cos Y 


2 uX 


- 2 tanY 


(B.5) 


dY 


1 


dy WCOSY u 
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and 


dp. 


w 


dy 

= - P„tan\ 

dp 

u 

2Pu>^ PV 

dy 

E*wcosY u^ 

dp 

^Y 

w p 2 p 

w u 


- 2np^) 


dy 


, ^ X sinV 

cos^Y cos^Y E^wcos^Y ^Y 


(B.6) 


respectively. Now, to maximize the Hamiltonian (B.4), we consider 
the switching function 


$ = E'" p^ - 2up^ (B. 7) 

Then, for to be maximum with respect to we use \ = \ if 

^ max 

$ > 0, and we use \ = \ . if ^ < 0, In the finite time interval during 

min ® 

which = 0, we have X. = variable. Since for maximum range glide, 
in the plot for \ in Figure 9 there is an interval in which the optimum 
X is variable and near unity, we have the approximate singular 
relation 

$ = E*p - 2up = 0 (B. 8) 

^Y n 

By taking its derivative, using Eqs. (B.5), (B. 6), and (B. 8) itself, 
we have 

wp =2p [”1 + ^ (1 + E ^ )1 (B , 9 ) 

w u L E ‘ ^ J 


As the linear control does not appear in this first derivative, take 
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the derivative again. This time the linear control appears and is 
solved to be 


= 


w cos Y 


u 


(B. 10) 


This is the approximate but explicit control law in Eq. (6. 31), 

Now, let us use \i as an approximate solution to linearize the 
Hamiltonian (B.l). Then 




X =\i 


8^ 

9X 


X =X, 


(X - Xi) 


w cos Y 


, ,, , /wco 

= wp tanY + p ( 


u 


u E‘'wcosY 


- 2 tan Y 


' u 


( 


2p 


w cos Y 




(B. 11) 


The corresponding linearized state and adjoint equations are 


and 


dp 


w 


dw 

dy 


= wtanY 


du w cos Y 
dy " 


u 


dY 


E’^u 

X 


E"" w cos Y 


2 ^ 


1 


dy 


dy w cos Y u 


/cosY 

p tanY - p — 
w '^u V E u 


u 


E* w^ cos Y ^ ^ 




w cos Y 


(B. 12) 


dp. 


u 


/ w cos Y 
“ ^uV E*u' 


dy 

dy 




1 


\ 

) ■ u^ 


wp 


E ■' w cos Y 


^w / w sin Y u sin Y 

^Y ^ \ E*u ^ E*wcos^Y cos^Y/ 


(B. 13) 


u sin Y 


cos 


Py X sin Y 
w cos^ Y 



106 


respectively. Again, by putting the switching function equal to zero, 
we have 

p^ = 2 w p^ cos Y (B. 14) 

By taking the derivative of Eq. (B. 14), it gives 


E w p 


= Pu [' 



2 w^ cos'* V 
E"u^ 


- ^ wcos^Y sinY + ~ sinYl (B. 15) 

u w J 


Then if we take the derivative of Eq. (B. 15), we can finally solve a 
new explicit control law which has been given in Eq. (6. 32). 
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APPENDIX C 


Derivation of Formula (8.26) 



Figure C. 1. (6^, (|>p and (0^ , cj)^) 


We have 4*^' and (9^ , and we want (9^, By considering 

the right spherical triangle M^AI, we immediately have 


tana = 


tan 0£ 

cos 4'.' 
1 


(C. 1) 


and 


tan at = sin 0 ' tan ib.' 

I 1 


(C,2) 


Since bj = cj?^' - ai , taking the tangent of hi and using (C. 2) in it, we 


have 
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tanbi = 


tan4>^ - sin 9^' tan4J.’ 
1 + sin0^* tan 4>£ tan 4^.' 


(C.3) 


Now, from the right spherical triangle M. A I, 


cos I = cos 9* sin 4'.’ 
f 1 


(C,4) 


On the other hand from the other right spherical triangle M^BI, 


cos I = tanb cot b^ 


(C.5) 


Solving for (tanb) from (C,4) and (C.5) and using (C.3), 


tanb = 


cos 0^* sin4^^* (tan 4>£ - si 
1 + sin9 ' • ^ 


sin 0 ' tan 4'.') 
f 1 


I 

1^' tancjj^' tan4j£ 


(C.6) 


As 9^ = a + b, again taking the tangent and using (C. 1) and (C. 6), we 
finally have 


tancj)' sin 4'.’ 

tan 9 . = ^ + tan 9 J cos 4’.' 

f cos 9' f 1 


(C.7) 


This is the first formula in (8.26). 

For the second formula, from the right spherical triangle 


M^BI, 


sin 4>^ = sinbi sin 1 


(C.8) 


and from the right spherical triangle AI, 


cos 4^/ = cos ai sini 


(C. 9) 


Hence , 
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sincj)^ 


sin bi cos 4*^' 


cos a .1 


(C. 10) 


Now, taking the sine on both sides of bj = - aj and then using 

(C. 2) , we have 


sin (f)^ = sin <|j^‘ cos 4^/ 


sinQ* cos sin 4 '.* 
f f 1 


(C, 11) 


This is the second formula in (8. 26). 
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